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Abstract 



We present further evidence for the bulk nature of the phase transition line in the fundamental- 
adjoint action plane of SU(3) lattice gauge theory. Computing the string tension and some 
glueball masses along the thermal phase transition line of finite temperature systems with 
N t = 4, which was found to join onto the bulk transition line at its endpoint, we find that 
the ratio \fa /T c remains approximately constant. However, the mass of the ++ glueball 
decreases as the endpoint of the bulk transition line is approached, while the other glueball 
masses appear unchanged. This is consistent with the notion that the bulk transition line 
ends in a critical endpoint, with the continuum limit there being a 4 theory with a diverging 
correlation length only in the ++ channel. We comment on the implications for the scaling 
behavior along the fundamental or Wilson axis. 
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1 Introduction 



The study of the phase diagram of fundamental-adjoint pure gauge systems, 

S = Pf D 1 " ^ReTrf/p] + (3 a £[1 - ^Tr^TrC/p], (1) 

in the early 80 's revealed a non-trivial phase structure with first order (bulk) transitions in 
the region of small (3f [[I], In particular a line of first order transitions points towards the 
fundamental axis and, after terminating, extends as a roughly straight line of bulk crossovers 
to the fundamental axis and beyond. This non-trivial phase structure, and in particular the 
critical endpoint, has been argued to be associated with the dip in the discrete /3-function of 
the theory with standard Wilson action, which occurs in the region where the bulk crossover 
line crosses the fundamental axis. 

In a couple of recent papers Gavai, Grady and Mathur and Mathur and Gavai ||] 
returned to investigating the behavior of pure gauge SU(2) theory in the fundamental- 
adjoint plane at finite, non-zero temperature. They raised doubts about the bulk nature of 
the phase transition and claimed that their results were consistent with the transitions, for 
lattices with temporal extent N t = 4, 6 and 8, to be of thermal, deconfining nature, displaced 
toward weak coupling with increasing N t . On the transition line for a fixed N t there is then 
a switch from second order behavior near the fundamental axis to first order behavior at 
larger adjoint coupling. In a Landau Ginzburg model of the effective action in terms of 
Polyakov lines, Mathur || reported that he could reproduce the claimed behavior seen in 
the numerical simulations. These results, should they be confirmed, are rather unsettling, 
since they contradict the usual universality picture of lattice gauge theory with a second 
order deconfinement transition for gauge group SU(2). 

Puzzled by the finding of Ref. |J, we studied the finite temperature behavior of pure gauge 
SU(3) theory in the fundamental-adjoint plane || [5]]. We obtained results in agreement with 
the usual universality picture: there is a first order bulk transition line ending at 

= (4.00(7), 2.06(8)). (2) 

The thermal deconfinement transition lines for fixed Nt (being of first order for gauge group 
SU(3)) in the fundamental adjoint plane are ordered such that the thermal transition for 
smaller N t occurs to the left, at smaller (3f, than that for a larger N t . In this order the 
thermal transition lines join on to the bulk transition line. The thermal transition line for 
N t = 4 joins the bulk transition line very close to the critical endpoint. This is shown in 
Figure 4 of Ref. 0, reproduced here as Figure [l|. 

To solidify this picture, which is in agreement with the usual universality scenario, we 
have continued the investigation studying zero-temperature observables. We computed the 
string tension and the masses of some glueballs, in particular the ++ glueball, along the 
thermal transition line for N t = 4. For universal continuum behavior y/a and the glueball 
masses should be constant along the thermal transition line for a fixed Nt, leading to constant 
ratios T c j sfo and m g j y/a. 
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Figure 1: The phase diagram together with the thermal deconfinement transition points for 
N t = 2, 4, 6 and 8 from Ref. f7j. The lower plot shows an enlargement of the region around 
the end point of the bulk transition. 



2 



Of course, at small N t , corresponding to a large lattice spacing a, we expect to see some 
deviations from this constant behavior. However, we find large deviations for m ++ / \fo as we 
approach the critical endpoint of the bulk transition along the Nt — 4 thermal transition line. 
The scalar glueball mass decreases dramatically, much more than what could be expected 
from simple scaling violations at a large lattice spacing. On the other hand, this is not really 
a surprise since at the critical endpoint at least one mass in the ++ -channel has to vanish. 
We elaborate on our findings in the next sections and then discuss the implications for the 
scaling behavior along the fundamental (or Wilson) axis. 

2 Observables and Analysis 

We have made simulations of the model with action (H) along the thermal transition line 
for Nt = 4, and continued along the bulk transition line, on a 12 4 lattice. Experience 
indicates this size to be sufficient to avoid significant finite size effects. The simulations were 
carried out with a 10-hit Metropolis algorithm tuned to an acceptance rate of about 50%. 
Observables were measured every 20 sweeps after thermalization. Close to the fundamental 
axis this resulted in essentially statistically independent measurements, while closer to the 
critical endpoint the autocorrelation time was significantly increased. We computed finite T 
approximants to the potential from time-like Wilson loops constructed with 'APE'-smeared 
spatial links |8| to increase the overlap to the ground state potential. On and off axis spatial 
paths were considered with distances R — n, y2n, \f?>n and \^En, with n — 1, 2, ...an 
integer. The string tension was then extracted from the usual fit 

V(R) =V -^ + l (g l (R) - i) + aR. (3) 

Here Gl(R) is the lattice Coulomb potential, included in the fit to take account of short 
distance lattice artefacts. Our fits are fully correlated x 2- fits with the correlations estimated 
by bootstrap, after binning to alleviate autocorrelation effects. The results of the best fits 
are listed in Table [|. 

We also computed glueball correlation functions in the ++ , 2 ++ and 1 H channel that 
can be built from simple plaquette operators. In an attempt to improve the signals we built 
these plaquettes not only from the original, but also from the smeared links, already used 
for the computation of the potential. Not surprisingly for computations done around the 
critical coupling for the N t = A thermal phase transition, we did not obtain a significant 
signal in the l +_ channel and only an effective mass from time distances t = 0/1 in the 
2 ++ channel - we had 500 measurements everywhere, except for (3 a = 2.0, near the critical 
endpoint where the number was increased, as given in Table |I[ In the ++ channel we got a 
signal at distance t — 1/2 at small (3 a and out to t = 3/4 at j3 a = 2.0. Our best results are 
also given in Table |TJ 

The quantities ^/a and m ++ are shown in Figure |^ plotted versus [3 a . As can be seen, 
\fa remains approximately constant along the thermal transition line for N t = 4 - the errors 
shown in the figure are statistical only; no estimate of the error from the uncertainty in the 
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Table 1: The results in the neighborhood of the N t = 4 thermal transition line. Comments: 
(a) ^/a at /3 Q = 0.0 comes from ||, m ++ from fllPf ; we did not take their best value, but 
rather the effective mass from the same distance as at f3 a = 0.5; rri2++ is from JTI]]; (b) in the 
disordered phase and (c) in the ordered phase on the bulk transition. 

determination of /3f c has been attempted except for f} a = 2.0. There, the computation has 
been repeated for two nearby couplings, also listed in Table [I]; the variation with Pf becomes 
so rapid that the error in the determination of /3f c becomes the dominating factor. 

While our estimate for m 2 ++ , albeit an unreliable estimate since we had to use distances 
t = and 1 to obtain enough of a signal to extract an effective mass, also remains approxi- 
mately constant, m ++ shows a remarkable decrease as the critical endpoint is approached. 
This observed behavior suggests that the mass in the ++ channel vanishes at the critical 
endpoint, thereby giving strong additional evidence for the existence of this critical endpoint, 
since at a critical point at least one mass, in the ++ channel to have a rotationally invariant 
continuum limit, must vanish. 

Since no other observable seems to be dramatically affected by the critical endpoint, 
we conjecture that the continuum theory one would obtain there is simply the (trivial) 4 
theory. To substantiate this claim somewhat, we made a fit to the scalar mass of the form 

m 0++ = A (P* a - p a f (4) 

expected near a critical point. A 3-parameter fit gave A = 0.76(11), p = 0.35(20) and 
P* = 2.02(6) with a x 2 of 0.29 for 2 dof. Note that the estimate for P* is in agreement with 
the previous estimate (0) obtained in from fits to the jump in the plaquette across the 
bulk transition line. Within its large error, the exponent p is compatible with the mean field 
value 0.5 of 4 theory, up to logarithmic corrections. Since the errors of the fit parameters 
are rather large we also made a fit with P* held fixed at its value 2.06 of (§). This fit gave 
A = 0.71(3) and p = 0.44(5) with a \ 2 of 0.47 for 3 dof. Again, p is compatible with the 
mean field value. Finally, a fit with P* = 2.06 and p = 0.5 both held fixed gave A = 0.68(1) 
with a x 2 °f 1-50 for 4 dof. This last, still very acceptable fit is included in Figure 0. 

While we believe to have assembled impressive further evidence for the bulk nature of 
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Figure 2: y/a (octagons) and m ++ (squares) as a function of (3 a along the thermal transition 
line for N t = 4. At f3 a = 2.25 we show the results from both phases at the bulk transition. 
The dashed vertical line gives the location of the critical endpoint, (|2|), with the dotted lines 
indicating the error band. The curve is a fit to m ++ = A(2.06 — Pa) 1 ^ 2 - 
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Figure 3: The potential in both phases on the bulk transition line at (f3f, f3 a ) 
on 12 4 lattices (□ and x) and on a 16 4 lattice (o). 



the phase transition line in the fundamental-adjoint plane and of its critical endpoint, we 
decided to make one further test. We computed the potential in both phases on the bulk 
transition line at f3 a = 2.25. This is about the location where the N t = 6 thermal transition 
line joins onto the bulk transition line according to Ref. 0. Hence, even in the ordered phase 
a lattice of size 12 should be large enough to obtain a reliable potential. To make sure that 
this is indeed the case, we repeated the computation in the ordered phase on a 16 4 lattice. As 
can be seen in Table p] the finite size effect on the 12 4 lattice does indeed appear negligible. 
On both sides of the bulk transition - at identical couplings: we did not observe even an 
attempt at tunneling from one phase to the other - we find a clearly confining potential. 
Therefore the bulk transition does not, as expected, affect confinement, provided the lattices 
are large enough. 

Of course, the string tension (and glueball masses) jump from one side of the bulk tran- 
sition to the other. Indeed, as the bulk transition line is approached from the fundamental 
axis the string tension varies more and more rapidly - and the thermal deconfinement transi- 
tion lines for different N t come closer together - i.e., the dip in the step /3-function becomes 
deeper until a jump is developed through the bulk transition line. 



3 Implications for Scaling 

In the previous section we have corroborated the existence of a first order bulk phase transi- 
tion line ending in a critical endpoint. We have provided evidence that physical observables 
are little affected by this phase transition line, except for a stronger dependence on /3f at 
fixed p a - a deepening of the dip in the step /3-function - and eventually a jump across the 
bulk transition line. The notable exception to this is the glueball mass in the ++ channel 
which decreases as the critical endpoint is approached and vanishes there. 

The influence of the critical endpoint on the ++ glueball appears still visible in the 
crossover region on the fundamental (Wilson) axis. This has first been argued in Ref. ||12|| . 
The ++ glueball is lighter in the crossover region, leading to a different scaling behavior 
than other observables. This can be seen in Figure [| where we show the latest data of 
Tj yfa from Ref. [|13|] and m ++ / yfa with ^Ja taken from Refs. || [L4j] and the glueball mass 



from Refs. JTD, [TJ| [T6|]. While T c /s/a stays approximately constant in the P interval shown, 
m ++ 1 \fa decreases visibly in the crossover region around P = 5.7. 

It has long been known that the scalar glueball mass scales differently in the crossover 
region than T c and the string tension, whose scaling behavior as a function of P deviates 
from asymptotic scaling but agrees with the step /3-function found in MCRG computations. 
The scalar glueball seemed more compatible with asymptotic scaling. However, in view of 
our findings that the different behavior of the scalar glueball mass comes from the influence 
of the nearby critical endpoint this asymptotic scaling behavior appears to be accidental. 
It would be interesting to determine the scaling behavior of other (pure gauge) observables. 
We conjecture that they would behave more like T c or ^/cr than like the scalar glueball mass. 
Unfortunately no reliable, large volume, glueball masses for J PC quantum numbers other 
than ++ for couplings in the crossover region exist in the literature to check this conjecture. 
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Figure 4: 10T C / yfa (octagons) and m ++ / ^fo (squares) as a function of (5 for the fundamental 
(Wilson) action. 
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